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(Nl ■ ABSTRACT 



' In an attempt to develop a deterministic theory for planet formation, we examine the 

accretion of cores of giant planets from planetesimals, gas accretion onto the cores, and 
their orbital migration. We adopt a working model for nascent protostellar disks with 
a wide variety of surface density distributions in order to explore the range of diversity 
[ among extra solar planetary systems. We evaluate the cores' mass- growth rate Mc 

through runaway planetesimal accretion and oligarchic growth. The accretion rate of 
cores is estimated with a two-body approximation. In the inner regions of disks, the 
^ ' cores' eccentricity is effectively damped by the ambient disk gas and their early growth 

■ is stalled by "isolation". In the outer regions, the cores' growth rate is much slower. If 

some cores can acquire more mass than a critical value of several Earth masses during the 
persistence of the disk gas, they would be able to rapidly accrete gas and evolve into gas 
giant planets. The gas accretion process is initially regulated by the Kelvin-Helmholtz 
contraction of the planets' gas envelope. Based on the assumption that the exponential 
decay of the disk-gas mass occurs on the time scales ~ 10^ — 10'' years and that the disk 
mass distribution is comparable to those inferred from the observations of circumstellar 
disks of T Tauri stars, we carry out simulations to predict the distributions of masses 
and semi major axes of extra solar planets. In disks as massive as the minimum-mass 
disk for the Solar system, gas giants can form only slightly outside the "ice boundary" 
at a few AU. But, cores can rapidly grow above the critical mass interior to the ice 
boundary in protostellar disks with 5 times more heavy elements than those of the 
minimum-mass disk. Thereafter, these massive cores accrete gas prior to its depletion 
and evolve into gas giants. Unimpeded dynamical accretion of gas is a runaway process 
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which is terminated when the residual gas is depleted either globally or locally in the 
form of a gap in the vicinity of their orbits. Since planets' masses grow rapidly from 
lOM® to lOOM®, the gas giant planets rarely form with asymptotic masses in this 
intermediate range. Our model predicts a paucity of extra solar planets with mass 
in the range lO-lOOM^ and semi major axis less than 3AU. We refer to this deficit 
as a "planet desert". We also examine the dynamical evolution of protoplanets by 
considering the effect of orbital migration of giant planets due to their tidal interactions 
with the gas disks, after they have opened up gaps in the disks. The effect of migration 
is to sharpen the boundaries and to enhance the contrast of the planet desert. It also 
clarifies the separation between the three populations of rocky, gas-giant, and ice-giant 
planets. Based on our results, we suggest that the planets' mass versus semi major axes 
diagram can provide strong constraints on the dominant formation processes of planets 
analogous to the implications of the color-magnitude diagram on the paths of stellar 
evolution. Finally, we show that the mass and semi major axis distributions generated 
in our simulations for the gas giants are consistent with those of the known extra solar 
planets. Future observations can determine the existence and the boundaries of the 
planet desert in this diagram, which can be used to extrapolate the ubiquity of rocky 
planets around nearby stars. 

Subject headings: extra solar planets - solar system: formation - stars:statics 

1. Introduction 

More than 100 extra solar planets have been discovered so far. This sample is sufficiently large 
that some clues and constraints on the formation of extra solar planets may be inferred from their 
kinematic properties. Figure 1 shows the distributions of semi major axes (a) and masses (Mp) of 
the known extra solar planets as of July 2003 (http://www.obspm.fr/encycl/encycl.html). Various 
radial-velocity surveys may have exhaustively found most of the planets in the solar neighborhood 
with radial velocity amplitude ^ 10 ms~^ and relatively small semi major axes, a < 4-5AU, 
which correspond to orbital periods < 10 years. The mass distribution in these regions appears to 
follow a weak power-law function in Mp with an abrupt upper cut-off at around 10 Mj, (Jorissen 
et al. 2001; Tabachnik & Tremaine 2002; Zucker & Mazeh 2002) where Mj is a Jupiter mass and 
is ~ 32OM0 (M® is an Earth mass). 

But, Figure 1 also shows a correlation between the Mp and a distributions. Zucker & Mazeh 
(2002), Patzold & Rauer (2002) and Udry et al. (2003) discussed a deficit of massive (Mp > 5Mj) 
and short period (a < 0.2AU) planets. Here we focus on another domain of paucity in the Mp — a 
distribution: the regions a > 0.2AU and Mp < Mj (Figure 1). This paucity is not due to the 
observational selection effects of precision and legacy baseline since planets in this domain have 
Vr > 10 ms~^ and orbital periods P less than a few years (Udry et al. 2003). In this paper, we 
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suggest that this paucity reflects the condition for and the process of planetary formation. 

We interpret the observational data on Jupiter-mass extra solar planets with the conventional 
core accretion model of gas giant planets (Wuchterl et al. 2000 and references therein). In the 
standard model, terrestrial planets and solid cores of gas giants are formed through the coagulation 
of planctcsimals with initial sizes ~ 1-10 km, e.g., (Safronov 1969; Wcthcrill 1980; Hayashi ct al. 
1985). Cores which are heavier than a few lunar masses attract the disk gas and form hydrostatic 
atmospheres. But, if a solid core reaches a critical mass of > several (§3.1), the pressure 
gradient in planet's gaseous atmosphere can no longer support against its gravity and the planet's 
atmosphere collapses onto its core (Mizuno 1980; Bodenheimer k, Pollack 1986; Pollack et al. 1996; 
Ikoma et al. 2000). Subsequent disk gas accretion onto the solid core forms a gas giant planet. 

The cores' asymptotic masses and their accretion time scales are determined by surface density 
of the planetesimal swarm and orbital radius (Kokubo & Ida 2002). In general, the formation of 
the gas giants is favored slightly exterior to the "ice boundary" at a few AU (Kokubo & Ida 2002), 
where the surface density of dust particles increases by a factor three to four as ice condenses. This 
enhancement of the dust surface density facilitates the formation of large cores. Well beyond the 
"ice boundary" (> lOAU), the cores' growth rate is so small that the residual gas in the disk tends 
to be depleted before it can be accreted. In very massive protoplanetary disks, however, cores can 
reach the critical mass interior to the ice boundary (Lin & Ida 1997; Kokubo &; Ida 2002). Since 
the core accretion is relatively fast in the inner regions of the disk, gas accretion onto cores can 
start well before the gas is depleted. Unimpeded dynamical accretion of gas is a runaway process 
which is terminated when the residual gas is depleted either globally or locally in the form of a gap 
in the vicinity of the gas giants' orbits. But, prior to a self-limiting growth stage, the planets grow 
so rapidly from lOM® to more than lOOM®, that they rarely emerge with an asymptotic mass in 
the intermediate range (10 — IOOMq). We will refer to such a planet-deficit domain in Mp versus 
a distribution as a "planet desert" . 

After their formation, gas giant planets with sufficiently large masses to open up gaps also 
undergo type II orbital migration which is well coupled to the viscous evolution of their nascent 
disks (Lin &; Papaloizou 1985). If these planets' tidal and magnetospheric interaction with their host 
stars can halt their inward migration, there would be an accumulation of "short-period planets" 
(equivalently "close-in planets", or "hot Jupiters"), near the disk inner edge (< O.IAU). Dozens 
of short-period planets have been discovered (Marcy et al. 2000) and their masses appear to be 
relatively small (Udry et al. 2003). These planets have either formed with relatively low masses 
(see §4) or lost some of their initial mass through Roche lobe overflow (Trilling et al. 1998; Gu et 
al. 2003). The latter process is only effective in the regions close to the stellar surface. Massive 
short-period planets may also be more vulnerable to post formation orbital decay due to their 
enhanced tidal interaction with their host stars (Zucker &; Mazeh 2002; Patzold k, Rauer 2002). 

Gas giant planets may attain a in the range of at 0.2 — 3AU, only if the disk gas is timely 
depleted during their migration (Trilling et al. 2002; Armitage et al. 2002). Type II migration 
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proceeds at a faster pace for planets with relatively lower mass and smaller orbital radius. As 
shown in §5, the intermediate-mass planets can open up gaps in inner regions (a < lAU) and easily 
migrate to the proximity of the host stars. Larger planets tend to migrate less rapidly than the 
intermediate-mass planets. 

Intermediate-mass planets with insufficient masses to open gaps also interact tidally with 
the disk gas through their Lindblad and corotation resonances (Goldreich &: Tremaine 1980). In 
principle, any imbalance in the torque can lead to "type I" migration (Ward 1986). The effect 
of the corotation resonances (Tanaka et al. 2002) and secondary instabilities in the vicinity of the 
coorbital region (Balmforth &; Korycansky 2001) may be able to quench the magnitude and sign 
of the net torque. Observationally, there is no evidence that type I migration has brought in any 
water-filled terrestrial planets into the inner Solar system. 

Well outside the ice boundary (> 3AU), the intermediate-mass planets can slowly accrete 
nearby planetcsinials. Prior to the depletion of the gas, however, intermediate-mass planets are 
unable to open up gaps in outer regions. In the absence of type I migration, they cannot migrate 
to the inner region from larger orbital radii. Cores in the outer regions of the disk can attain 
masses in excess of IMq after the gas depletion (§2). The mass distribution for long-period planets 
are expected to be continuous. But, with the residual gas depleted, the intermediate-mass planets 
which formed lately cannot migrate to fill the planet desert. 

Based on these considerations, we suggest that the intermediate-mass planets may be rare and 
their deficit gives rise to a "planet desert" at 0.2 — 3 AU. In §2, we present the model of core 
accretion and predict cores' masses as functions of their nascent disks' mass and semi major axis. 
In §3, we include the effects of gas accretion onto the cores. Assuming an appropriate distribution 
of disk masses that may be comparable to the observationally inferred mass distribution, we carry 
out, in §4, Monte Carlo calculations to simulate the mass and semi major axis distributions of 
extra solar planets. The predicted distributions clearly show a paucity of planets with masses 
~ 10 — lOOAf® at a < 3AU. §5 discusses the effects of radial migration of giant planets due to their 
tidal interaction with a gas disk. We show that type II migration sharpens the boundaries and 
enhances the contrast of the planet desert. In §6, we summarize our results. 

2. Core Accretion 

In this section, we briefly recapitulate the current theories for solid core accretion. Our main 
objective is to determine the cores' growth rate and asymptotic masses. 
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2.1. Core accretion model: the feeding zone 



In an early stage of planetesimal accretion, the largest planetesimals have the shortest mass- 
doubling time scale. This process is commonly referred to as the 'runaway growth', e.g., (Greenberg 
et al. 1978; Wetherill & Stewart 1989; Aarseth et al. 1993; Kokubo & Ida 1996). In the limit that 
their velocity dispersion is relatively small, field planetesimals are readily captured by nearby 
runaway bodies (cores). However, the differential Kcplerian speed of the planetesimals prevents 
them from reaching the distant runaway cores until they have acquired sufficient orbital eccentricity. 
The region within which that the runaway cores can directly accrete planetesimals is commonly 
referred to as the 'feeding zone'. When the runaway cores capture all the planetesimals in their 
feeding zones, their growth stops as they acquire an 'isolation mass'. Lissauer (1987) evaluated the 
isolation mass as a function of the initial surface density of the solid components in protostellar 
disks. Assuming a negligible velocity dispersion for the planetesimals, the full width of the feeding 
zone is given by Aoc ^ 2^/12t'h — 7rH, where 



is the Hill radius (Roche lobe radius) for a solid core with a mass Mc around a host star with a 
mass M*. 

However, when a core becomes sufficiently massive to excite the eccentricity and to pump up 
velocity dispersion of the neighboring planetesimals, its feeding zone expands. The increases in 
the relative speed between the planetesimals and the cores also reduce the effective cross section 
of the cores such that their growth slows down (Ida & Makino 1993; Aarseth et al. 1993; Rafikov 
2003). The relatively massive cores are more effective in stirring up the velocity dispersion of their 
neighbors so that the reduction of their effective cross section is larger. Consequently, the growth 
time scale becomes an increasing function of the cores' mass. Such a self-regulated growth results 
in the formation of a population of comparable-mass cores which are embedded in a swarm of 
small and relatively eccentric planetesimals (Kokubo &: Ida 1998, 2000; Thommes et al. 2002). In 
this system, the orbital spacings of the cores (equivalently, full width of their feeding zone, Acc) 
are nearly equal to one another and are as large as ~ lOrn- This maximization of spacing is the 
combined consequence of the distant perturbations between the cores and the effect of dynamical 
friction from the small planetesimals on the cores (Kokubo & Ida 1995). Kokubo & Ida (1998) 
called this coagulation process, the 'oligarchic growth'. Because the width of the feeding zone is 



during the oligarchic growth stage, the isolation masses are slightly larger than those estimated by 
Lissauer (1987). 

The concept of an isolation mass induced by the depletion of the feeding zone is based on the 
assumption that the orbits of the cores do not evolve with time. In reality, the cores undergo radial 
diffusion due to their scattering of the field planetesimals, e.g., (Murray et al. 1998; Ida et al. 2000) 




(1) 



Aoc ~ lOrn 



(2) 
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and their tidal interaction with the gas disk (type I migration), e.g., (Goldreich & Tremaine 1980; 

Ward 1986). Under the action of gas drag, the small planetesimals and dust particles may also 
undergo migration, e.g., (Adachi et al. 1976). These effects can lead to an expansion of the feeding 
zones. Thus, in a prescription for the asymptotic mass of cores' growth, it is worthwhile to consider 
the possibility that the cores' feeding zone may be extended throughout the disk. 



2.2. The cores' growth rate 

The growth rate of the cores is determined by the velocity dispersion (a) or equivalently the 
eccentricity distribution of the planetesimals and cores. Eccentricity damping due to their tidal 
interaction with the disk gas (Artymowicz 1993; Ward 1993) inhibits the cores from attaining any 
significant growth in their velocity dispersion. In contrast, the eccentricity of the low-mass planetes- 
imals is less effectively damped. Instead, it is excited by cores through recoil of dynamical friction. 
Through a series of 3D numerical simulations, Kokubo & Ida (2002) showed the planetesimals' 
eccentricity is larger than the ratio of the cores' rn and their semi major axes so that the cores' 
growth time scale is well described by a simple two-body approximation. 

In the modest to high-cr limit, the cores' accretion rate at an orbital radius a is (Safronov 1969) 



where Or = (GM*/a=^)V2 

is a Kepler frequency, pd and Ed are the spatial and surface density 
of solid components (pd ~ EdSlK/c), and Rc is a physical radius of the core. In order for the 
planetesimals within Aoc to reach the cores, 

— ~ Aoc ~ 10 — ^ a (4) 



so that the solid cores' accretion time scale is 

where Tk = 27r/f^K [= (M,/M0)-V2(a/lAU)3/2 years] is the Kepler period. By taking into account 
the effect of gas drag on the planetesimals' velocity dispersion, Kokubo k, Ida (2002) derived a more 
detailed expression (see their eqs.l5 and 16), 

Tcacc- 1.2X10 l^^Qg^^_2j (^lAuJ \M^) \Mq) 



llAuJ 



2.4 X 103gcm-2y VIAU/ \m^^g 



(6) 
years. 
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where Sg is the surface density of the gas in the disk, m is the mass of typical planetesimals, and 
Mq is the Sun's mass. Since rc,acc = Mc/Mc, we find from eq. (6) 



Mcit) 



VlAuJ V 



TT- Ma 

Mq 



3.5 X lO^years / \10gcm~2y \^2.4 x lO^gcm^^y yiO^^g 



(7) 



where we neglect variations of and Sg with time. This approximation is adequately satisfied in 
the inner regions of the disk where Tc^acc 

is small compared with the gas depletion time scale. 

In the outer regions of the disk, the rate of mass growth is reduced from the above expression 
during the depletion of the disk gas. In eq. (48), we introduce a prescription for the global depletion 
of the gas. For this prescription of Sg, the evolution of the core mass becomes 



\3/v^ \3/ \6/5/ \ -2/5 



X 



1.4 X lO^yearsy V10gcm"V V2-4 x 103gcm-2y VlO^^S/ 



where Tdisk is the disk depletion time scale. In the limit of t <C Tdisk, eq. (8) reduces to eq. (7). 

According to eq. (8), the cores attain an asymptotic mass. But for t > raisk, the assumptions 
for circular orbits of the cores and damping of planetesimals' a by gas drag break down. When all 
the gas is depleted and all the heavy elemental content is retained by the oligarchic cores, their a 
increases to become comparable to their characteristic surface speed a ~ T^urf = {GMd RcY^'^ and 
their collisional cross section reduces to their physical cross section. In this limit, the growth time 
scale for the final planetary embryos become 

(9) 

where Mc, Rq, and pc are a mass, radius, and an internal density of the embryo. Here, we refer to 
solid planets after depletion of the gas as "embryos" rather than "cores". The reduced collisional 
cross section significantly lengthens the growth time scale for the embryos from that of the cores. 
The corresponding mass for the embryos is 

M,(t, .M.(.,^,+ (^19S±-)' (^)' (^)" {^y^M,. (10) 



2.3. Disk structure models: the temperature and surface density distributions. 



The cores and embryos' mass, as expressed in both eqs. (8) and (10), is a function of the 
surface density of the heavy elements in the form of condensed grains and planetesimals, Ed- The 
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distribution of is poorly determined theoretically and observationally. In principle, the 

distribution is determined by the infall pattern of the heavy elements, the condensation, growth, 
and sublimation rates of grains, and the particles' interaction with ambient disk gas. Here, we 
adopt an ad hoc power-law distribution such that 



based on a phenomenological minimum-mass solar nebula model (Hayashi 1981). The power-law 
dependence may be the consequence of self-similar evolution. We use this prescription as a fiducial 
model. We plot as a function of a for = 0.1, 1, and 10 in Fig. 2. Using these models, we assess 
the formation properties of planetary systems as a function of S^. The step- function variable r/icc is 
introduced to express the effect of ice condensation/sublimation across the ice boundary (aice ^ a 
few AU; see eq. [15]). For the inner disk with a < aice, ?7ice = 1 and for outer regions with a > aice, 
Vice ~ 3 — 4. Following the minimum-mass solar nebula model, we adopt Tjic^ — 4.2 to represent the 
Solar abundance distribution (Hayashi 1981). 

We also introduce a scaling parameter for the total disk mass, /d, which is assumed to be a 
constant parameter throughout the disk. In the case of the minimum-mass disk model for the Solar 
system (Hayashi 1981), the scaling parameter ~ 0.7. This value would be larger if the formation 
of the solid cores and their retention into the Solar system planets were not highly cfRcicnt. Radio 
observations of the continuum radiation in the mm-wave length range measure the dust emission 
from the disks. The observationally inferred total mass of dust, M^, in the protostellar disks around 
classical T Tauri stars ranges from IO'^Mq to 3 x lO^^M©, e.g., (Beckwith Sz Sargent 1996). The 
disk mass determinations are somewhat uncertain due to the poorly known radiative properties of 
the grains. Nevertheless, their differential dust content provides a reasonable evidence for a greater 
than an orders of magnitude dispersion in Mj. In addition, the image of the disks is not well 
resolved in many cases. A rough magnitude of /d can be inferred from the the total mass of the 
disk under the assumption that all disks have similar sizes (a few tens AU to a hundred AU). The 
observationally inferred mass of T Tauri disks corresponds to a range of /d ~ 0.1 — 30. We note 
that disks with /d substantially larger than unity do commonly exist. 

The presence of gas affects the velocity dispersion and therefore the growth rate of the cores. 
It is also important for the formation of gaseous giant planets. In principle, the magnitude of the 
gas surface density Eg is determined by the gas infall rate and it evolves as a consequence of the 
viscous stress in the disk, e.g., (Lin &: Papaloizou 1985). Similar to all astrophysical accretion disks, 
the effective viscosity is generally assumed to be due to turbulence. But, in regions of protostellar 
disks where planets emerge, the cause of turbulence is highly uncertain (Papaloizou Sz Lin 1995). 
For computational simplicity, we also assume that Eg has a similar power law dependence on a as 
Ed such that 



We vary the scaling parameter /g to obtain different disk models and to test the dependence of the 
planetary properties on the mass of the gas disk. This assumption implies that the relative ratio 




(11) 




(12) 
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of gas (volatiles) to dust (refractories composed mainly of heavy elements) 

(13) 



Eg _ /g240 



/d Vice 

is uniform throughout the disk except across the discontinuity at the ice boundary. While /d//g = 
1 would represent the Solar abundance, a different choice of /d//g would represent the relative 
enrichment /depletion of the heavy elements in the disks. 

In the above equations, we introduce a separate mass scaling parameter /g for the gas in order 
to consider the possibility that the solid components may evolve independently from the gas. The 
basic concept of the minimum-mass solar nebula model is based on the assumption that all the 
heavy elements are retained while the original disk gas was depleted. Theories suggest that gas 
diffuses under the action of viscous stress, but dust and planetesimals migrate under the influence 
of the gas drag effect such that they tend to evolve independently from each other. Gas may also 
be preferentially depleted through the photo evaporation e.g., (Shu et al. 1993; Matsuyama et al. 
2003) and winds while dusty grains may be preferentially retained by the disk (Shang et al. 2000). 

At the moment, we can observe the dust component around most, and molecular hydrogen 
around a few, young stars. Among the classical T Tauri stars, there is no apparent dependence of 
the mass of the dust disks on the host stars' ages. For stars with ages greater than ~ 10® — 10^ 
years, the dust mass inferred from the mm-wave length continuum is much reduced (Beckwith & 
Sargent 1996; Wyatt et al. 2003). The near-IR signatures of disks also evolves on a similar time 
scale (Haisch et al. 2001). This decline has been interpreted as an evidence for dust growth and 
planetesimal formation rather than the depletion of heavy elements (D'Alessio et al. 2001). In this 
paper, we adopt the conjecture that the surface density distribution of the dust plus planetesimals 
does not change except in those regions where they have been totally accreted by the cores. 

Very little information is available on the evolution of Eg and /g, but there is no indication 
of divergent depletion pattern between molecular hydrogen and mm-size dust emission (Thi et al. 
2001). In addition, the presence of spectral evidences and UV veiling for ongoing accretion are 
well correlated with the mm flux for the classical T Tauri stars. Both signatures of protostellar 
disks also disappear together among weak line T Tauri stars (Duvert et al. 2000). Based on these 
observations, we adopt the most simple-minded assumption that gas is depleted on a similar time 
scale even though the gas and dust depletion proceed through different mechanisms. 

In eq. (11) we introduce a parameter T^ice to take into account of the possibility of ice conden- 
sation. In order to determine the radius of the ice condensation, Oice, we need to specify the gas 
temperature in the disk. In most planet-forming regions of interest, the gas density is sufficiently 
large to maintain a thermal coupling between the gas and the dust particles. In the optically thin 
regions of the disk, the dust and the gas are heated to an equilibrium temperature (Hayashi 1981) 
such that 

T = 280 i^r'" (h.) K = 280 ""^ (Ml.) K. (14) 
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For solar-type stars, the stellar luminosity L* is roughly proportional to Mf. Since a-jce corresponds 
to the radius where T ~ 170K, we find that 

/ M \ ^ 

Hereafter, we adopt Oice = 2.7AU (L* = Lq) until §4.1. In the Monte Carlo calculations in §4.2 
and 5, we consider variation of Oice- Since scale height of a gas disk is /i ~ Cg/^K where Cg is sound 
velocity, from eq. (14), we find 




2.4. The cores' asymptotic mass 
Based on the distributions of given by eq. (11), we find from eq. (7), 

by assuming m = lO^^g. Equivalently, the core accretion time scale is (eq. [6]) 

Note that the actual time scale to reach Mc is 3rc,acc but not Tc^acc- The numerical factor of 3 

1/3 ' ' 

comes from the dependence that Tc acc 

The unimpeded core growth which is described by eq.(17) can not be sustained indefinitely. 
In the limit that damping effect due to disk-planet interaction (Artymowicz 1993; Ward 1993) 

regulates a small velocity dispersion for the cores, their growth terminates when they consume all 
the planetesimals in their feeding zones. If the radial diffusion/migration of the cores is neglected, 
the isolation mass (Lissauer 1987; Kokubo k Ida 1998, 2000) 

Mc,iso = 27raAacSd (19) 

represents the asymptotic mass of the cores prior to the depletion of the disk gas. If we adopt 
Attc = lOrn (where Mc in rn is replaced by Mc^so); the isolation mass would be 



^c4so=^0.16 — — - — ^ — ^ Me. (20) 



10gcm-2y VlAUy \10ru) V^e 

Upon acquiring Mc^iso, the cores consume all the nearby planetesimals. These cores cannot 
cross each other's orbits, until the gaseous-disk mass has decreased to < lO^^Me, equivalently 
/g < 10~^, and the damping effect due to disk-planet interaction becomes weak enough (Iwasaki 
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et al. 2002; Kominami &: Ida 2002). Even in a gas free environment, a swarm of planetesimals has 
a tendency to evolve into a quasi-static oligarchic state. Kokubo & Ida (2002) performed three- 
dimensional N-body (with N=10,000) simulations of planetesimal accretion in gas-free disks with 
a wide variety of surface density distribution of the planetesimal swarm. Their results indicate 
that the process of planetary accretion is fully consistent with the oligarchic growth scenario, i. e. 
the cores initially emerge with a characteristic orbital separation of Aoc ~ lOrn- In later stages, 
however, Aoc increases to ~ 15rH, which may be due to the radial diffusion of the cores (they 
neglected radial migration due to disk-planet interaction and did not include the effect of gas drag 
on a population of small planetesimals). Here we adopt Auc = lOrn- Prom eq. (11), we estimate 
the core masses in protoplanetary systems with eq. (20) such that 

M,,i,o O.lGrii'fJ' iTXuf ^^^^ 

If the effect of radial diffusion/migration of the cores is included, the final core mass could be- 
come larger than Mc^iso- In the limiting case that a migrating core can acquire all the planetesimals 
interior to its orbital semi major axis, it would attain a maximum asymptotic mass 

M,..„ ~ = 1.2 (^^) (^)' M,. (22) 

Based on the distributions of given by eq. (11), we find 

Mcnoiso - l-2??ice/d (^) M®. (23) 

Growth beyond isolation is also possible after the gas depletion and the elimination of its 
damping effects (Iwasaki ct al. 2002; Kominami & Ida 2002). After the gas depletion, the velocity 
dispersion of the residual cores grows until they cross each other's orbit. The final stages of cores' 
growth may proceed through giant impacts (Lissauer &; Stewart 1993). Eventually a few surviving 
embryos acquire most of the residual planetesimals and less massive cores during the late oligarchic- 
growth stage. Since their collisional cross section is reduced to their geometrical surface, the mass 
of the embryos evolves with the growth rate given by eq. (9), 

dMe Me 



r .4xl0VT'r^fr^V'7^V'Y^V'' years ' ' 

The asymptotic embryos' masses are given by eq. (20) with Aoc ~ Vsmi/^K, 

M,,^.0.&^{^y''[^)\,. (25) 



Such a scenario is consistent with the current model for the origin of the Moon (Hartmann 
& Davis 1975; Cameron &; Ward 1976), and the apparently young isotropic age of the Earth, in 
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comparison with asteroid parent bodies (Halliday et al. 2000). Note that even if large-mass cores 
can be assembled through giant impacts, they can no longer initiate gas accretion because the disk 
gas needs to be severely depleted prior to their emergence. The limited gas supply during the 
late phase of core growth through giant impacts is also consistent with the intermediate mass and 
mostly-ice composition of Uranus and Neptune. 

In outer regions, the growth of the cores may be terminated before their mass reaches -Me^iso 
or Mc,noiso- When their characteristic surface speed V^urf = ^JGMe/Re exceeds their local escape 
speed from the host star, I4scp = \/2GM^Ja, the embryos' 90" grazing scatterings of nearby bodies 
can lead to the bodies' ejection (Thommes et al. 2003). In this high a limit, the ratio of collision to 

ejection probabilities is /cap ~ (Kjscp/Kiurf)*^- When /cap become significantly less than unity, the 
growth of the embryos is terminated. Based on previous numerical simulations (Lin &: Ida 1997), 
we set /cap ~ 0.1 such that most scatterings lead to ejection and the embryos' asymptotic mass 
becomes 

-^e,sca 2/cap 

a ^^^^ 

where is the embryos' physical radius. For an average ice embryo, with an internal density 
Pe ~ Ig cm~^, the asymptotic mass of the embryos is 

This limiting mass at a 20 — 30AU is comparable to those of the cores of Uranus and Neptune. 
Note that this limit only applies to the scattering of planetesimals and cores and it is not a limiting 
mass for gas accretion (see §3). For notations of masses, time scales, etc., see Table 1. 

In principle, the largest cores may emerge with sufficient asymptotic masses to initiate gas 
accretion (see below). But in many regions, this large mass can only be attained after the gas is 
depleted. Thus, they are prohibited from evolving into gaseous giant planets. In this paper, we 
consider the formation of the gas giants only if some cores are able to acquire the critical mass for 
the onset of rapid gas accretion before the disk gas is depleted. We assume that the coalescence 
between the cores, through giant impacts after the gas depletion, do not lead to substantial further 
growth in the mass of their gaseous envelopes despite the possibility of substantial increases in their 
core mass. 



2.5. Core growth at different orbital radii 

With the above prescriptions, we estimate, with eq. (7), the cores' growth rate and asymptotic 
mass limit. The initial mass of the cores is arbitrarily set to be a small value, = lO^'^g. The 
choice of an initial value of does not affect the results, since t ^ acc increases with Mq . We plot , 
in Figure 3a, the core masses, after 10^ years and 10^ years, as a function of at various locations 
a. For this determination, we assume that the growth of the cores is truncated with a mass Mc^iso 
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(eq. [21]). We also assume Sg/S(j = 240. A similar plot with the cores' truncation mass set by 
Mc^noiso (eq. [22]) instead of Mc^iso is shown in Fig. 3b. 

As shown in eq. (7), the cores' growth is relatively fast in regions where is high and a is 
small. The results in Figures 3 show that for this region of the parameter space, the cores' growth 
is terminated with their asymptotic masses within 10^ years. In contrast, the growth is relatively 
slow for regions of parameter space with relatively small and large a so that Mc{t) remains to be 
very low in the right lower regions in Figures 3. An "accretion wave" propagates from the top-left 
to the lower-right regions as time passes. In the very large Sd and large a regions, Aoc ~ lOrn 
for Mc^iso exceeds the magnitude of a. Thus, the appropriate expression for the truncation mass is 
noiso rather than -Mc in these regions. 

Next, we consider planetary systems with a distribution similar to that derived from a 

minimum-mass solar nebula model (see eq. [11]). Figures 4 show some examples of the time 
evolution of the cores' masses for (a) /d = 1 and (b) = 10. The jumps in the cores' masses at 
2.7AU in these figures reflect a transition in the values of S^. These results indicate that in the 
inner regions of the disk, the gas depletion time scale is longer than the build-up time scale for the 
cores and all the field planetesimals are accreted onto the isolated cores. But, for the minimum- 
mass solar nebula model (with /d ~ 1), iWc^igo at 0.7 and 1 AU is considerably smaller than the 
present masses of Venus and Earth (see eq. [21]). These planets could have acquired their present 
mass through either giant impacts after depletion of the gas or the accretion of planetesimals prior 
to the gas depletion if they have diffused or migrated over large radial distances. In this case, the 
cores' asymptotic mass would be Me,iso (eq. [25]) or Mc^noiso (eq. [23]) which, at 0.7 and lAU, is 
comparable to the mass of the Earth and Venus. 

For a more comprehensive parameter study, we show in Figures 5, the dependence of the cores' 
mass on a and /d after 10^ years and 10^ years. In these calculations, /g does not decline and is set 
to be always equal to /d since the stages prior to depletion of the gas are considered. The magnitude 
of Sd for the minimum-mass solar nebula model corresponds to /d — 1. On the gas depletion time 
scale ~ 10^ - lO'^ years, the most massive cores attain a critical mass (~ several M^) just outside 
the ice boundary (2.7AU). Although the magnitude of both -Mc,noiso and Mc,iso increases with a, 
the expression in eq. (17) indicates that much beyond aice, the cores do not have sufficient time to 
attain their asymptotic values prior to the depletion of the disk gas. In this case, the formation of 
gas giants within a minimum-mass solar nebula is most likely to be initiated just exterior to the 
ice boundary. But around many classical T Tauri stars, the disk mass inferred from the mm-wave 
length continuum observation corresponds /d ^ 5 (see §1). The results in Figs. 5a and b suggest 
that in very massive protoplanetary disks with /d > 5, cores can rapidly attain masses greater than 
the critical mass for gas accretion (see next section) even inside the ice boundary. 

After the gas depletion, Mc^iso expands to Mc^noiso or iWe^iso- Afe,sca can also limit embryo 
growth (see, §4.2.3 and Fig. 10). In the calculations in the following sections, we set the asymptotic 
masses of cores and embryos as follows: 
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• the isolated core model 

- mm(Mc,iso, Afc,noiso), when /g > 10~^ (before the gas depletion). 

- min(Mc,iso,Mc ,noiso; 

Me 

,sca)) when /g < 10 (after the gas depletion). 

• the non isolated core model 

- Mc,noiso, when /g > 10~^. 

- min(Mc,noiso,Me,sca) when /g < 10~^. 



Up to now, we only consider the drag effect induced by the gas on the damping of the plan- 
etesimals' and cores' velocity dispersion. In this section, we consider the accretion of gas onto 
sufficiently massive cores. The necessary condition for any gas accumulation near the cores is that 
their surface escape velocity must exceed the sound speed of the disk gas. In regions of the disk 
where dust can condense to form planets, the sound speed is less than a few km s~^ such that even 
lunar-mass cores can attract gas to their proximity. However, the disk gas first accumulate into 
a static envelope around these cores. In this envelope, the gravity of the cores is balanced by a 
pressure gradient which is maintained by the energy release during the planetesimals bombardment 
onto the cores. 



But, when the cores' mass (Mc) becomes larger than a critical value (Mc^crit)) the state of 
quasi hydrostatic equilibrium can no longer be maintained and a phase of gas accretion onto the 
core is initiated (Mizuno 1980; Bodenheimer Sz Pollack 1986). Thereafter, the planetesimals and 
gas accretion proceed concurrently, so that the planetary growth can no longer be expressed in 
an analytical form as that in eq. (17). Based on the existing numerical results we adopt a simple 
prescription to evaluate Mc^crit and the gas accretion rate, dMp^g/dt. 

The critical core mass depends on the cores' rate of planetesimal accretion (Mc) and the opacity 
(k) associated with the disk gas. Faster accretion and higher opacity (relatively large Mc and k) 
result in a warmer planetary atmosphere and the enhanced pressure gradient stalls the gas accretion 
(Stevenson 1982; Ikoma et al. 2000). Based on a scries of numerical models, Ikoma et al. (2000) 
find that the critical core mass has the following dependence. 



3. 



Gas accretion 



3.1. 



Gas accretion onto cores 




(28) 
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We adopt here a power index 1/4 for the dependence on Mc- If the dust particles have a similar 
distribution and abundance as the interstellar medium, k ~ Icm^g^^. The actual magnitude of k 
is uncertain since the dust particles' size distribution and the fraction of heavy element in the gas 
and solid phases are poorly known (Podolak 2003). For computational simplicity, we neglect, the 
dependence of -Mc^crit on k in our prescription below. 

The gaseous envelope of proto giant planets contracts on a Kelvin Helmholtz time scale when 
their masses Mp (including both their solid cores and gaseous envelope) become larger than Mc,crit- 
Ikoma et al. (2000) derived through numerical calculation that 

'''' 

with power index 6 ~ 8 and c ~ 2.5. However, with different opacity table, the values of h and c 
can be changed to be ~ 10 and ~ 3.5 respectively (Ikoma, private communication). Bryden et al. 
(2000b) obtained 6 ~ 10, c ~ 3.0 by fitting the result of Pollack et al. (1996). We adopt here 6 = 9 
and c = 3.0 and also neglect the dependence on k in our prescription below. 

The contraction of the gas envelope reduces the local pressure gradient. Consequently, the gas 
from the surrounding disk flows toward the cores to replenish the envelope at a rate 

dM^^M^ (30) 

dt TKH ^ ' 

This rate is well below the Bondi accretion rate (Bondi 1952) when Mp ~ Mc^crit so the gas accretion 
rate is insensitive to the boundary condition in the disk. The magnitude of dMp^g/dt due to the 
Kelvin Helmholtz contraction increases rapidly with Mp. In principle, dMp^g/dt is limited by the 
Bondi accretion rate (Bondi 1952), the Keplerian shear in the gas inflow (Tanigawa &; Watanabe 
2000), mass diffusion rate through the disk Bryden et al. (1999), and the availability of gas near 
the orbit of the accreting planet. But, these effects become important only when Mp has already 
become a significant fraction of Jupiter's mass Mj. For Mp ~ Mj, tkh is reduced to the dynamical 
free-fall time scale and the growth time scale Mp/(dMp^g/dt) for Bondi accretion is reduced to 10^ 
years. Thus, the total gas accretion time for the protoplanets to attain their asymptotic masses is 
essentially determined by tkh at Mp ~ Mc,crit- 

Note that when Mc reaches Mc^iso -^c is reduced to zero. Consequently, Mc^crit also vanishes 
and gas accretion can start for any value of M^. However, if 

/ \ -1/3 

Mc<Me,KH = Mp(rKH = rdisk)-4.6f^^^J M®, (31) 

the gas accretion rate is so small (eq. [30] ) that Mp is not significantly increased by the gas accretion 
during the first 10^ years after they were formed. In order for the gas giants to actually form and 
acquire Mp Mc^crit, M: must be larger than Mc^kh in addition to the requirement of Mc > Mc^crit- 
Both of conditions can be satisfied for Mc > 4M0, if the cores' rate of planetesimal accretion is 
significantly reduced from that in eq. (3). 
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3.2. A prescription for the growth of gas giant planets 

Based on the above considerations, we introduce the foUowing simple prescription for the 
purpose of numericahy integrating the formation and early evolution of the gaseous giant planets. 
These approximation captures the essence of the planet formation process. 



Proto planetary growth due to planetesimal accretion is computed with the following rate 
equation, 

= = M^{M^ = Mp), (32) 



dt 

where Tc,acc is given by eq. (18). The value of Mc is obtained from eq. (3) with the core mass 
Mc being replaced by the protoplanet's total mass Mp which also includes the mass of the 
gaseous envelope. 

During the growth of the cores, we check whether Mp is less than the critical core mass for the 
onset of gas accretion. Based on the discussions in the previous subsection, the magnitude of 
M;,crit is approximated with a simplified version of eq. (28) such that 

In principle, as the field planetesimals arc being accreted onto the cores, their surface density 
in the feeding zone declines, so that both Mc and Mc,crit decrease. To take this effect into 
account, we use time-dependent values of Mc calculated in (32) and set Mc,crit = for cores 

with Mc > Mc,iso or Mc,noiso- 

When Mc exceeds Mc^crit, gas accretion onto a core is set to start. The gas accretion rate is 
regulated by the efficiency of heat transfer and it is approximated by a simplified version of 
eqs. (29) and (30) such that 



dt TKH 

M, 



-3 (34) 



The planets' mass Mp includes both their gaseous envelope and cores. Due to the expansion 
of the growing planets' feed zones, planetesimal accretion onto the cores continues after the 
gas accretion has started. But, we neglect the collisions between cores until after the disk gas 
is depleted as we have discussed in the last section. 

The above gas accretion rate does not explicitly depend on the ambient conditions. The 
approximation is appropriate provided there is an adequate supply of the disk gas. However, 
the disk gas may be depleted either globally through processes such as the viscous evolution. 
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outfiow, and photo evaporation or locally through gap formation. These effects are taken 
into account in the determination of the asymptotic mass of the giant planets (see the next 
subsection) but not in their accretion rates. 

• In the limit that the growth of the gas giant planets is terminated by gap formation, we also 
consider their dynamical evolution with or without the type-II orbital migration (see the next 
section). 

3.3. Termination of gas accretion 

Equation (34) shows that the gas accretion rate rapidly increases with Mp. However, there are 
several asymptotic limits to the unimpeded runaway gas accretion process. These limits arise due 
to the local and global depletion of disk gas. 

3.3.1. Local gas depletion 

In an inviscid disk, gas accretion onto the cores is terminated when they consume all the gas 
in their feeding zone and attain an isolation mass. The isolation mass for the gas giants can be 
obtained by replacing in eq. (20) with Eg such as 

where the Hill (Roche lobe) radius rn is evaluated with Mp. For a fiducial magnitude of Mg^jgo, we 
used the width of gas feeding zone Aog = 2rH to be a nominal value. Since it flows normally on 
circular orbits, gas which is much beyond Aog cannot reach the accreting cores without modifying 
their specific angular momentum. Using the Sg distribution in eq. (12), we find 

M...Wr(^)-(^)^'^(^)-'%., 

which is comparable to Mj at the present location of Jupiter. Since Sg and equivalently /g decrease 
with time, Mg also decreases with time. dM^^g/dt is set to be zero when Mp becomes larger than 

Mg,iso. 

However, the observed accretion flow from protostellar disks to their host stars (Hartmann 
et al. 1998) requires mass and angular momentum transfer throughout the disks. In a viscously 
evolving disk, gas can continuously diffuse into the feeding zone until the growing planet attains a 
mass (Lin & Papaloizou 1985), 

M.,.. . ^M. . 40„ (^) \4 . 30 (^) (^) (^) (37) 



-18- 



where we used eq. (16) and a-prescription for the effective viscosity v (Shakura & Sunyaev 1973), 

in which u = ah'^Q^ where a is a dimensionless parameter. For Mp > Mg^visj an embedded planet 
induces a strong tidal torque to open a gap near its orbit, provided its Hill (Roche lobe) radius is 
larger than a critical value rH,c = h (Lin &: Papaloizou 1993; Bryden et al. 2000a). We find that 
''H = ''H.c when 

, / a \ 3/4 o / a \ 3/4 / \ 

M„,=.3.8xlO-(_) M. 1.2 X 10^ (_) (38) 

Because Mg th is generally larger than Mg^vis except for some very large a cases, gas inflow onto 
the planet would actually be halted, when Mp reaches Mg^th- In some previous hydrodynamic 
simulations, embedded protoplanets continue to accrete after the above condition is satisfied albeit 
at a reduced rate (Artymowicz & Lubow 1996). In oder to take this possibility into account, we 
also modify, in one set of calculations, the truncation condition to rH,c = 1-5^. 



3.3.2. Global gas depletion 

Before their Mp increases to the isolation mass, Mgas,iso, protoplanets can accrete at the rate 
given in eq. (34). But, for Mp > Mgas.iso, additional disk material can no longer be brought into 
the feeding zone by the expansion of its boundary. Nevertheless, in a viscously evolving disk, the 
feeding zone may be replenished. We also consider the extreme case that the protoplanets' tidal 
interaction does not impinge the flow of the gas from the disk to their feeding zone. 

Although gas diffusion can replenish the feeding zone, the rate of gas accretion onto the pro- 
toplanet is limited to 

<'^P.-n,inf^,M,.,) (39) 



dt \TKH 

where the mass transfer rate in the disk, 

Mdisk = 27rSgaFdisk, (40) 
the gas diffusion speed is determined by the efficiency of viscous transport 

With the a prescription for a steady disk, Vdisk ~ —ah'^ClK/a. 

The magnitude of a in the planet formation region is uncertain. We assume that the mass 
transfer rate throughout the disk is comparable to an observational inferred accretion rate from 
protostellar disks onto classical T Tauri stars (Hartmann et al. 1998) which can be fitted by 



Mdisk ~ Mo(t/Tdisk)"^-^ 



(42) 



where Mq = IO^^Mq year"-*^ for Tdisk = 10^ yr (Calvet et al. 2000). The actual data have an 
order of magnitude spread over this fit. The accretion rate onto the planets is limited by the gas 
replenishment rate of their feeding zone, Mdiskj when their 

^ > Mdisk (43) 

TKH 

or equivalently when their Mp becomes larger than 



t 



-3/8 



(see eqs. [29] and [43]). The limiting replenishment rate only modifies dMp^g^s/dt for protoplanets 
more massive than both Mg^gup and Mg^igo. Prom eqs. (36) and (49), we find that these conditions are 
satisfied at the advanced stages of disk evolution when the surface density of the disk is significantly 
depleted. 

At any given location, the emergence of cores which can dynamically accrete gas occurs at an 
epoch 

1 /3 

Tonset = 3re,a.c(Me = M,,ent) ~ 8 X W'rjr^f-'/' (^) {^^) y^^'^ (^5) 

for M* = IMq. According to eq. (42), the amount of mass can be supplied into the feeding zone 

for t > Tonset IS 



Mg,feed = Mdisk(t)c/t ^ 1.5Mdisk(i = Tonset)Tdisk - 60 (^)'^' /d'^'° (loXu) ^'^^^^ 



Tdisk 



10''' years 
(46) 

If such supply is considered, the isolation mass can increase from Mg^iso to Mg^iso+Mg^feed- However, 
at large a, Mg^feed is not a large enough. Since sufficiently massive cores (with Mc > Mc^crit) take 
longer to form in the outer regions of the disk, the mass transfer rate can limit the amount of 
material supplied to the feeding zone, Mg^feed, to less than the mass of Jupiter. 

The gas accretion is ultimately limited by the amount of residual gas in the entire disk. For 
our disk model, the maximum available mass is 

When Mg^noiso < -Wg,th) Mp is truncated by Mg^noiso rather than Mg th- This condition occurs when 
/g < 3(a/lAU)^/^(M*/M®)^/^. Similar to Mg^iso, Mg^noiso exphcitly depends on the magnitude of 
/g, i.e., the amount of residual gas in the disk. 

As we have already indicated in §2 that the gas content are observed to decline concurrently 
with the dust (but not the planetesimal) content in protostellar disks. Since the near IR signature 
for dust in the inner region declines together with the mm indicators of dust in the outer regions of 
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protostellar disks (Duvert et al. 2000), we assume that the residual gas in the these disks is depleted 
globally in a self similar manner. Wc introdTice a prescription to approximate the evolution of Sg 
during the depletion phase of the disk such that 



where Sg^o is an initial distribution (eq. [12]) and Tdisk is gas depletion time scale. Equivalently, 
decreases with time from its initial value /g,o as 



Based on the observed properties of protostellar disks (Beckwith &; Sargent 1996; Wyatt et al. 
2003), we set Tdisk = 10^ years. 

In order to limit the extent of the parameter studies, we also assume here that /g,o = fd, that 
is, the gas-dust ratio is initially equal to the Solar composition. We denote both /g,o and /d to be 
/disk- We consider the cases with /disk = 0.1 — 30 as mentioned in the last section. We do not vary 
the value of /disk during the cores' growth but it is locally set to zero when the cores have reach 
their asymptotic masses. 

Based on these prescriptions, we set an asymptotic mass for the gas giants. We assume that 
the gas accretion onto it is terminated when a protogiant planets either consumes all the nearby 
residual gas in the disk or has acquired a sufficiently large mass to open up a gap near its orbit. 



In this section, we consider dynamical properties such as the mass-semi major axis distribution 
which planets are formed with. This distribution evolves with the orbital migration of the gas giant 
planets. We will discuss effect of orbital migration in the next section. 



We compute a series of models with Tdisk = 10^ years. We determine the asymptotic (at 
t = 10^ years) magnitude of Mp, including that of the gaseous envelope if any, as a function 
of /disk (= fd = /g,o) Ei'iid o (see Figures 6). For these models, we adopted the isolated core 
model (§2.5). We adopt five different prescriptions for the planets' asymptotic mass after their gas 
accretion is truncated: (a) Mg^iso, (b) Mg^noniso, (c) Mgas,vis, (d) Mg^th, and (e) Mg^iso + Mg^feed- 
For the evaluation of the asymptotic mass of the gas giants, we adopt Aog = 2rH in (a) and (e), 
M* = IMq in (c) and (d), and a = 10^^ in (c). Since the result of (e) is almost the same as that 
of (d), it is omitted. There are some uncertainties in the amount of gas diffusion to the planet's 
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4. Nascent Orbital Distribution of Potoplanets 



4.1. Conditions for the emergence of gaseous and ice giant planets 
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feeding zone. These five upper mass limits provide a range in the asymptotic mass of the planets. 
In Figs. 7, we show the results of a second series of models computed with the non isolated core 
model (§2.5). All other parameters are identical to the first series of models. 

For disks with a uniform /disk ~ 1 — 3, there are three distinct regions in Figures 6. Similar 
regions also exist in Figures 7 for disk with uniform /disk ~ 0.3 — 1. In the inner (a < a[cc) regions, 
the asymptotic mass of the cores Mc is smaller than 4M^, which is the minimum mass (Mc^kh) 
needed to undergo a significant amount of gas accretion within 10^ years. Since gas is depleted on 
this time scale, gas giants can not form in these region. Although after gas depletion, the cores' 
mass can exceed several through giant impacts between them, they can no longer accrete gas. 
These massive cores evolve into terrestrial planets similar to the Earth and the mass distinction 
between the terrestrial planets and gas giants is preserved. 

In the outer regions of the disk where a > 10 — 20AU, the cores' planetesimal accretion rate 
is very low. Although the asymptotic mass limit for the cores is large compared with several M^, 
there is inadequate amount of time for the cores to attain a significant amount of mass prior to 
the depletion of the disk gas. Nevertheless, low mass cores can slowly accrete gas because, for 
very low planetesimal accretion rates Mc, the values of Mc^crit is also small. But, in the low-Mc 
limit, the gas disk is severely depleted by the time Mc becomes comparable to Mc^kh- In advanced 
stages of gas depletion, the magnitude of both Mg^iso and Mg noiso becomes smaller than Mc and 
the cores cannot acquire a massive gaseous envelope. This limited supply of residual disk gas may 
have quenched the ability of Uranus and Neptune to attain a massive envelope (< 2M®) despite 
their present ice-core mass being > ISM® (Pollack et al. 1996). We refer to these heavy element 
enhanced, late comers with modest gaseous envelopes as ice giant planets. 

In the intermediate region slightly exterior to the ice boundaries, an abrupt increase in the 
cores' asymptotic mass allows the cores to initiate the onset of rapid gas accretion. Gas giants with 
Mp > IOOMq are preferentially formed in this bright region in Figures 6 and 7. The boundaries 
of the region where gas giants can form are set by the requirements that: (a) Mc > 4M0 and (b) 
3Tc^acc ^ Tdisk- The results in Figures 6 and 7 indicate that conditions (a) and (b) are satisfied, for 
/disk ^ 1, at atg (the "rocky/gas-giant planets' boundary") > 3AU and Ogi (the "gas/ice giants' 
boundary" ) < 10 — 20AU. Terrestrial planets form interior to otg and ice giants form exterior to 
Cgi. The values of of atg and Ogi are expressed in eqs. (53) and (54) below. 

We now consider atg and agi for more general values of /disk- Transition from cores to gas 
giants requires gas accretion and can only occur prior to gas depletion. But, in the presence of 
gas, the core growth is limited by their asymptotic masses. If the asymptotic mass is determined 
by the isolation process, it would become Mc^iso- In this case, condition (a) can only be satisfied if 
Mc,iso > 4M0. From eq. (21), this necessary condition corresponds to 
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Similarly, in the non isolated limit, the condition that Mc^noiso > corresponds to 

/disk > ftg - ^rd (iXu) ^ (^1) 

(see eq. [23]). These estimates agree with the terrestrial/gas-giant planets' boundary in Figures 6 
and 7, respectively. In comparison with the results in the first series, the emergence of the more 
massive cores in the second series of simulations enhances the probability and speeds up the onset 
of gas giant formation for modest values of a. 

In condition (b) (Sr^acc ^ "^disk); the quantity 3rc,acc is the actual growth time scale for the 
cores (see eq. [18]). Prom eq. (18), the condition for 3rc,acc — "^disk corresponds to 

/disk > /gi - 4 (— ) (^^J [W^s) (lOAuJ • 

This condition does not depend on whether the asymptotic mass of the cores is limited by -M^ jgo 
or Mc,noiso- For r/ice = 4, Mp = 4M^, and Tdisk = 10^ years, this equation indicates that /gi 2± 
4(a/10AU)V2, which also agrees well with the gas/ice giants' boundary in Figures 6 and 7. 

As the analytical derivation shows, the two boundaries separating the three regions nurture 
the formation of terrestrial planets, gas and ice giants and they do not depend on conditions of 
truncation of gas accretion. In the outer regions of the disk, the gas/ice giants' boundary is also 
insensitive to the asymptotic mass of the cores. At large a, the barrier for reaching Mc^crit is 
imposed by the persistence of disk gas on time scales Tdisk rather than the asymptotic limit of the 
core's mass. 

For disks with arbitrary values of /disk; the boundaries 

_ / 64/d^sk^icc AU (the isolated core case) 
I 9/disk^ice (the non isolated core case). 



»»f(5^)-°'"f^V'°'7 -.^V'°'^'aU (54) 



and 

„ . _ c;f 14/27 fl 

4 J \4:Mq)J VlO^ears^ 
are set by the conditions /disk = ftg and /gi respectively. In massive disks with considerable heavy 
elemental contents, /disk may be considerably larger than unity. In these massive disks, Otg is less 
than the ice boundary. In that case, it is possible to form gas giants interior to the ice boundary. 
In disks with /disk smaller than a critical value 

/disk,min = 0.9or0.4 — — — [TTyr—-] , (55) 



4 / \4:M(sJ \W years J 

flgi < o-tg (the factors 0.9 and 0.4 correspond to the isolated and non isolated cases, respectively). 
These results are in agreement with the numerical results presented in Figures 6 and 7. 
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Since cores must attain at least 4:M^ before they can rapidly accrete gas, in metal deficit 
disks with /disk < /disk.min; the gas giants cannot form. For a minimum mass solar nebula model, 
/disk = 1 and is < /tg inside the ice line and > /gi outside the ice line such that gas giants can 
only form at large radii. S ince /disk,min is Smaller beyond the ice line, it sets a more stringent 
condition with rjice = 4. They are also consistent with the apparent paucity of metal deficit stars 
with planets, e.g., (Gonzalez 1997; Santos et al. 2001; Murray &; Chaboyer 2002) and (D. Fischer, 
private communication). 

In the gas giant formation domain, the total mass of the planets increases rapidly through 
runaway gas accretion after it exceeds Mc^kh ~ 4M0. Gas accretion continues until the gas is 
depleted either globally or locally near the planets' orbits. Since the planets' growth time scale 
from the onset of rapid gas accretion to their asymptotic mass is relatively short compared with the 
disk depletion time scale, planets with Mp ~ 10 — lOOM® occupy a very limited domain interior to 
a < 3AU, as shown in Figures 6 and 7 (with the exception of the case where Mp is limited by Mg^yis)- 
In outer regions (a > 3 — lOAU), however, ice giants in such mass ranges can be formed. Their 
cores attain sufficient mass to accrete gas only after the gas has already been depleted severely. 



Using the same theoretical model of core growth and gas accretion, we performed Monte Carlo 
calculations to produce a theoretical prediction for the Mp — a distribution of extra solar planets to 
compare it with observations (Figure 1). Planets are assumed to form with equal probability per 
interval of log a. We assume the /disk distribution as in Fig. 8 to be consistent with observational 
data (Beckwith &; Sargent 1996; Wyatt et al. 2003). This is a gaussian distribution in terms of 
logio /disk with a center at log^o /disk = 0.25 and dispersion of 1. 

We omit the high /disk tail at /disk > 30, since such heavy disks are self gravitationally unstable 
in outer regions. The gravitational stability parameter (Toomre 1964) for our disk models is 



The mass transfer associated with the gravitational instability would rapidly decrease local surface 
density Eg to a value such that Q > 1 (Lin & Pringle 1990; Nakamoto & Nakagawa 1994). Surface 
density in outer regions (a > lOAU) in massive disks (/disk ^ 20) may be adjusted. In principle, we 
should use the adjusted surface density distribution. However, since most of giant planets found 
in our model are formed at a < lOAU, we adopt a simple cut-off at /disk = 30. We also examined 
a more artificial distribution in which /disk is uniform in log scale from 0.1 to 30. But, the results 
are almost the same. 

Tdisk is distributed from 10^ to 10^ years uniformly in log scale, which is also consistent with 
observational data (Beckwith &; Sargent 1996; Wyatt et al. 2003). The range of stellar mass is also 
considered, so that Oice is changed as eq. (15). The location of Oice is essential for the formation 



4.2. Prediction for mass and semi major axis distribution 




(56) 
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of gas giants, as shown before. In order to directly compare with the observed data, we adopt a 

similar selection criterion for in the range of 0.7 — 1.4A^q. Correspondinglyj dice — 

1.3-5.3AU. 

We also assume a uniform distribution in log scale for M*. In each run, 10,000 initial conditions 
are taken. 

We carry out calculations for planetary formation in disks with an evolving gas content. In 
these models, the cores' accretion is truncated when their masses have attained asymptotic values of 
-^c,iso (the results with Mc^noiso are similar). These results are shown in Figure 9. Three truncation 
asymptotic masses resulting from the truncation of gas accretion are considered: (a) Mg^iso, (b) 
and (c) Mgth, and (d) Mg^yis- We adopt Aog = 2rH in (a), the critical Hill's radius rH,c being h 
and 1.5/i in (b) and (c), respectively, and a = 10~^ in (d). 

In Figure 9, the green filled circles and blue crosses represent rocky and icy planets with an 
insignificant amount of gaseous envelopes around them. The green and blue open circles represent 
gas-rich rocky and icy planets with gaseous envelopes which are one to ten times more massive 
than their cores. The red filled circles represent gas giants with gaseous envelopes more than ten 
times more massive than their cores. The result in panel (a) yields most of gas giants form with 
masses an order of magnitude larger than Jupiter. Even larger masses would be obtainable if the 
gas accretion is unimpeded until Mp becomes Mg^noiso or Mg^iso + Mg^feed- This inconsistency with 
the observed mass distribution of the extra solar planets is removed by the thermal gap formation 
condition as represented in panels (b) and (c). However, the viscous gap formation condition in (d) 
places an over strict upper limit on Mp so that gas giants cannot form. Numerical simulations of 
disk-planet interaction indicate that Mg^th is the appropriate upper limit for masses of gas giants 
(Bryden et al. 1999). 



4-2.1. Planetary desert 

In panels (a), (b), and (d), these results show a deficit of planets with intermediate mass 
(IO-IOOM0) and semi major axis less than 3 AU. The ranges of the deficit is different in panels 
(a), (b), and (c). The lower mass boundary of the planet-deficit domain in Figure 9 is demarcated 
by Mp ~ for all models. This value of Mp simply corresponds to the minimum core mass 

that can satisfy both Mc > Mc^crit and > Mc,kh (§3.1), which are the condition for the cores to 
undergo rapid gas accretion. Transition across this boundary is weakened by the growth through 
giant impacts after the disk depletion. The truncation conditions of the cores and gas do not affect 
this boundary. Future observational determination of this lower mass boundary will provide useful 
information and constraints on the critical core mass which can lead to the onset of rapid gas 
accretion. 

The upper mass boundary of the planet-deficit domain in Figure 9 is demarcated by the 
truncation condition for gas accretion. First, we consider the case of Mg^iso. The upper mass 
boundary of the planet-deficit domain corresponds to the minimum asymptotic mass of gas giants. 
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As shown in eq. (50), in the case of core accretion being truncated by Mc^iso, rapid gas accretion 
would occur if /disk ^ ftg — 8(a/lAU)^^/^. In the inviscid-gaseous disk model, the gas accretion 
process is terminated by isolation when the planets mass reaches a value which is given by eq. (36). 
Since fg^ = /disk in all the calculations described here, we obtain, by substituting /tg into eq. (36), 
a minimum truncated mass of gas giants in these regions to be 

^ 1.1 X lO^Me. (57) 

A similar calculation with /tg ~ 3(a/lAU)~-'^/^ (for cores to be truncated by Mc^noiso) yields 

~ 2.5 X lO^Me. (58) 

In regions of the disk where /disk < /tg; cores cannot undergo rapid gas accretion. But, in 
regions where /disk > ftg, cores accrete gas until Mp reaches M*^^^ x (/disk//tg)^^^- The process 
of gas accretion can be stalled by the depletion of the disk gas. But the time scale for the planets 
to attain their asymptotic masses is shorter than the gas depletion time scale so that only small 
fraction of planets may attain some intermediate masses. Therefore, most planets that undergo gas 
accretion have Mp > -'^^gasiso' ^^giso corresponds to the upper mass boundary of the planet- 
deficit domain. Note that the a-dependence vanishes in M*^^^, so that the upper mass boundary is 
essentially independent of a in panel (a). If asymptotic mass is regulated by Mg^iso, the observed 
upper boundary of the planet-deficit domain would provide constraints on the feeding zone width 
(Aflg) of gas giants and the critical surface density of solid components that can produce cores 
large enough for gas accretion (equivalently, /tg). 



4-2.2. Asymptotic mass of the gas giants 

Next, the results with gas truncation by Mg th are considered. These conditions are particularly 
appropriate since the ongoing accretion from protostellar disks onto nearly all classical T Tauri stars 
clearly reflects the consequence of viscous transport of angular momentum and mass in the disks. 
The thermal truncation condition does not include an explicit /disk dependence in the asymptotic 
mass of the gas giants. The mass dispersion of giant planets at given a is caused by the dispersion 
of the stellar mass M* and is relatively small. The giant planet mass as a function of a clearly 
reflects the gap opening condition (38). The upper mass boundary of the planet-deficit domain, or 
equivalently, the lower mass boundary of the giant planets domain, is demarcated by Mg^th with 
M* ~ O.TMq in the present models. The truncation condition for the core growth stage does 
not affect the results. Note that we have adopted some idealized truncation conditions. In many 
numerical simulations, e.g., (Artymowicz & Lubow 1996), there are indications that the gas flow 
on to the giant planets continues, albeit at some reduced rate, after the gap formation. Turbulence 
may be particularly important in sustaining the planets' growth until they have acquired larger 
masses than those shown in panel (b) of Figure 9. In order to demonstrate this dependency, we 
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consider in panel (c) an additional model with rH,c = 1.5/i. This relaxed condition yields 3 times 
larger Mp for the gas giants. 

The main difference between panels (b) and (c) is abundance of the green and blue open circles 
which denote the gas-rich terrestrial (rocky) planets and ice giants with the mass of the envelope 
do not exceed ten times that of the core (i.e. Mp < lOMc). Prom eqs. (21) and (38), we find 

M,,„ 0.75 X loVi'f ,r„"^ i^f (^) M.-- (59) 

Although rjice = 1 interior to the ice line, the critical value of /disk(— 17) for Mg^th < lOMc^iso in 
panel b (where rH,c = h) is within the range of /disk (0.1-30) we have adopted in our simulations. 
Outside the ice boundary where r^icc = 4, it is much more likely that for Mg^th < lOMc^iso- Even 
though the mass of the gas-rich terrestrial planets and ice giants are dominated by their gaseous 
envelopes, the metallicity of their envelope may be greater than that of the disk. But, in panel c 
where rH,c = 1-5^, Mg^th is much larger than Mc,iso- The asymptotic composition of the gas giants 
in these disks is diluted by the massive gaseous envelope. 



4-2.3. Radial extent of planet forming region 

The present scenario also indicates an outer boundary for the domain of gas or ice giant planet 
formation. The outer boundary for gas giants is set by Cgi (eq. [54]) where 3rc,acc(M:,KH) = Tdisk- 
Although in eq. (54), Ogj increases with /disk) it is limited to 30 AU even for the maximum value of 
/disk(= 30) we have adopted. This limited extent of the gas-giant formation zone only applies to 
the core-accretion scenario. In the gravitational instability scenario, gas giants are preferentially 
formed far from their host stars. The spatial distribution of gas giants may provide a discriminating 
test for the avenue of planet formation. 

After the gas depletion, the residual rocky and icy embryos can grow through giant impacts. In 
the above discussions, we have already indicated that core grow is limited by a) the availability of 
building-block material (through Mc^noiso; eq. [23]) and b) the retention efficiency (through Me^scai 
eq. [27]). But, even with unimpeded growth and infinite supply of planetcsimals, c) the maximum 
mass can be attained within the life span of the host star (r*), through collisional coagulation, is 
(eq. [101) 

M.,„. . Me(r., = 0.8 X 10^ [^J {^f (^)' {^Y' M,. (60) 
The boundary between conditions a) and b) is demarcated by noiso = -^e,sca it occurs at 

=3(511)-"' (^)"'\u. (61) 
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The boundary between conditions b) and c) is demarcated by Me,sca = -^c,iife ^^nd it occurs at 




Mc^lfc with t = IG years, Mc^noisoj and Mc^sca are plotted as functions of a for /disk = 1)3, 10, 
and 30. in Fig. 10. Rocky and icy planets can exit in the range surrounded by these masses. 
These analytic approximation is consistent with the upper limits of Mp shown in Fig. 9. In all 
panels of Fig. 9, the formation of both gas and ice giants is quenched beyond ~ 30 — 40 AU despite 
the availability of planetesimal building blocks. In panels a), b), and c), the gas-poor ice giants 
(represented by the blue crosses) do not extend above ~ 2OM0. Their growth is limited mostly 
by Mc^iifc. After the gas depletion, the enhanced velocity dispersion makes the accretion rate very 
low. In outer regions beyond ~ 30 — 40AU, this effect quenches ice giants' further growth. If the 
asymptotic mass of the giant planets is limited by the viscous condition (panel d), gas accretion 
beyond 10 AU would be quenched when Mp ~ 40 — lOOM®. Since their growth is speeded up by 
the gas accretion, they may be able to attain masses larger than Mg^yig over some intermediate 
range of a (3-10 AU) through coUisional coagulation after the gas depletion. 

For modest values of /disk such as that in the minimum mass nebula model (~ 1), atg < Cgi < 
dab < cbhc- In this case, rocky, gas and ice giant planets naturally form in order of increasing 
distance from their host stars. This expectation is consistent with the distribution of the planets 
in the Solar System. Outside Oabi the massive ice giants become the scattering agents which 
induce the residual planetesimals to be ejected to form cometary clouds and freely floating planets. 
This extrapolation is also consistent with the observed population of scattered Kuiper Belt objects 
(Duncan Sz Levison 1997). Beyond abc, the growth time is the limiting factor for the emergence of 
any sizeable protoplanets. In order to account for the absence of modest-size Kuiper Belt objects 
beyond ~ 40 — 50AU, /disk ^ a few is required. 



4.3. Comparison with observational data 

The present observational data (Figures 1) arc obtained with some selection effects. For 
example, only planets with radial velocity amplitude Vr > 10ms~^ and semi major axes a < 4— 5AU 
have been detected so far. The planets' mass inferred from these observations represent a set of 
minimum values and their real magnitude is a function of the inclination angle of their orbital plane 
to the line of sight. Many planets are found with a ~ 0.1 — lAU and Mp > Mj. Although with 
the model in which the gas accretion is truncated by first isolation, we can construct such massive 
close-in planets (Figure 9a), but their emergence requires massive disks which may be prone to 
gravitational instability. The second and third models (Figure 9b, c) do not produce Jovian-mass 
planets interior to ~ lAU because the disks' aspect ratio is relatively small. There are many 
possibility, within the framework of the core accretion scenario, to account for this discrepancy. 
Post-gap-opening accretion can lead to in situ formation of such planets. Orbital migration can 
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also bring in the relatively massive planets formed in the outer regions of the disks' aspect ratio is 
relatively large (see next section). 

The present data show an apparent lack of planets with Mp < 2OOM0 and a ~ 0.2 — 3 AU 
(see Fig. 1). When corrected for the inclination of their orbits, the demarcation boundary may 
be shifted to higher masses, by a factor of two, from that indicated in Figure 1. If this hint of a 
planet-deficit domain is confirmed by the future observations, it would provide useful constraints on 
the truncation condition for gas accretion. At their face value, the present data are consistent with 
the simulated results with a combination of the isolation condition with Mgas,iso and the thermal 
condition with Mg th- The viscous condition with Mg^yis limits planet masses to values much smaller 
than observed masses unless 1) an extremely large value of a is taken (which can be ruled out by 
the persistent time scale of protostellar disks) or 2) accretion continues after gap formation. 

Our numerical simulations also show an outer boundary for the planet-deficit domain at around 
3 AU. This boundary reflects the changes in the dominant process which limits the formation of 
gas giants. At a < 3AU, the condition of > is more important, because cores grow rapidly 
and reach their asymptotic mass prior to the era of gas depletion. But, at a > 3AU, the cores' 
growth is so slow, that it is more important for rc,acc ^ ^disk to be satisfied. For most cores that 
can start gas accretion, Tc,acc is of oder of Tdisk- Hence, many cores may not be able to initiate the 
process of gas accretion even though their mass can become larger than lOM® after gas depletion. 
As a result, planetary mass distribution at a > 3 AU is likely to be continuous and does not show 
a clear deficit. Thus, the outer a boundary of the planet-deficit domain indicates the critical semi 
major axis where core accretion timescales and disk gas depletion timescales are comparable. 

5. Orbital migration 

The existing data of known extra solar planets include a population of short-period planets 
with period down to 3 days. We consider the possibility of orbital migration in this section. 

5.1. Model of orbital migration 

There are several suggestions that, due to a torque imbalance, embedded protoplanetary cores 
can undergo rapid type I migration , e.g., (Goldreich & Tremaine 1980; Ward 1986). At 1 AU, the 
migration timescale for Mc ~ IM^ is only 10^ — lO^years (Ward 1997; Tanaka et al. 2002; Bate 
et al. 2003). If type I migration is considered, all the protoplanetary cores have a tendency to 
rapidly migrate to the proximity of their host stars, prior to gas accretion. Thus, type I migration 
appears to be inconsistent with prolific formation of extra solar giant planets as well as the efficient 
preservation of Earth mass planets in the habitable zones around the Sun. 

In contrast, type II migration naturally occurs when a planet acquires an adequate mass to 
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open up a gap. Gap formation process itself requires the dissipation of tidally induced density 
waves. Turbulent viscosity is a possible mechanism to provide such a dissipation mechanism. It 
can also regulate the rate of gas accretion onto the classical T Tauri stars. Once a protoplanet 
has acquired an adequate mass to open up a gap, its orbital evolution is coupled with the viscous 
evolution of the disk (Lin Sz Papaloizou 1985). 

In the next set of models, we include the effect of type II migration for gas giants which have 
already attained their asymptotic masses. For computational convenience, we introduce a simple 
prescription to follow the orbital evolution of the gas giant planets. The direction of the migration 
can be both inward and outward. We arbitrarily set an stopping semi major axis for the inwardly 
migrating planets to be ~ 0.04AU. There are several potential mechanisms for halting the migration 
near the stellar surface (Lin et al. 1996; Trilling et al. 1998), but none of these processes are effective 
beyond ~ 0.1 — 0.2AU. The migration can be halted beyond ~ 0.1 — 0.2AU by gas depletion (Trilling 
et al. 2002; Armitage et al. 2002; Matsuyama et al. 2003). However, both the planet formation and 
migration timescales at 0.2-3 AU is shorter than raisk ~ 10^ — 10^ years. We consider below the 
necessary condition to build up gas giants and deplete the disk during the course of their migration 
so that they would be relocated to this intermediate location. 

Here we consider only the type II migration in which orbital evolution of a planet embedded 
in a gap is coupled with the evolution of the disk (Lin et al. 1996; Trilling et al. 2002; Armitage et 
al. 2002). In a slowly evolving (quasi steady) disk, the angular momentum flux due to the viscous 
stress is approximately independent of the disk radius. For a simple approximation, we estimate 
the net angular momentum transfer rate, in the absence of any planets, at the radius of maximum 
viscous couple Rm, e.g., (Lynden-Bell &; Pringle 1974), 

Jm ^ ^SgZ/f^K,m-Rm - Y^g'"^m^K,m^m> (63) 

where the subscript "m" denotes values of various quantities at i?ni- 

Across the gap, the planets' tidal torque regulates the angular momentum transport. When 
the disk interior to their orbits is accreted onto their host stars or that exterior to their orbits 
diffuses to much larger radii, the embedded planets maintain the angular momentum transport 
balance (l/2)MpQK,papap ~ Jm by adjusting their orbital radius Op (Lin & Papaloizou 1985; Lin 
et al. 1996). This process leads to planetary migration on a rate, 

~ 3sign(ap - Rm)a^^ — ^ — J^K.m, (64) 

Op Mp \Ik,p V Op / 

where Ok,p is Kepler frequency of the planet's orbit at Op. 

The sign implies that after it has opened up a gap, an embedded planet at radius Op would 
migrate with the disk gas toward its host star if Op < Rm and away from the star if Op > Rm- 
Using eqs. (16) and (12) to substitute for Eg and h, we obtain a migration timescale, 

. ^ . 0.3 . ,„V- {^) (f ) (^)- y^. (e.) 
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5.2. An a-prescription for gas disk evolution 

The value of R^a depends on the Sg distribution. For the minimum-mass disk model, we 
estimate it to be ~ lOAU (Lin &; Papaloizou 1985). However, the observationally inferred size and 
mass of the disks have a large dispersion. A similar dispersion are expected in R^- During the 
viscous diffusion process, evolves outward if there is no infall to replenish the disk mass. In the 
absence of infall and outflow, Sg declines in most part of the disk as mass is accreted onto the host 
stars. From the conservation of the total disk mass 

Jtot ~ 7rSgi?^0K,m, (66) 

we expect 

k - -ft- <»^> 

If the disk depletion is mostly driven by photo evaporation or stellar wind ablation, the expansion 
rate of i?m would be slower. However, as shown below, in the model with Sg oc a^'^/-^ and h oc a~^/^ 
(equivalently, T oc a~^/^), planetary migration rate is independent of i?ni except for outer regions 
(a ^ Rm) where the direction of migration is outward. Hence, we only need an approximate 
prescription for the time evolution of Rm in our simulation. Corresponding to eqs. (48) and (67), 
we set 

Rrn = lOexp ( 1 AU. (68) 

VSrdisk/ 

In a quasi steady state, the radial velocity of the gas which is 

in regions well inside Rm- The scale hp is evaluated at Op. Note that |dp| cannot be faster than 
\ur\. Equivalently, Tmig cannot be smaller than 

For small planets embedded in relatively massive disks, Tmig is limited by T(jisk,acc- 

Also, when Op 

decreases to 0.04AU, migration is halted as mentioned above. 

Integrating dp with eqs. (65) and (68), we show examples of type II migration with our model 
in Figure 11, where time is scaled by Tmig at lAU (rmig,i). At a ~ lAU, planets can migrate 

1/2 

by a significant distance in the limit that rjep ^ "^mig,!- Since Tmig oc Op , the inner planets 
tend to migrate over a greater radial extent and the inward migration accelerates with time. In 
all figures, we indicate the evolution of Rm with dashed lines. Planets exterior to Rm migrated 
outward. Some planets which formed slightly outside the initial Rm migrate outward first. But, as 
Rm expands beyond their orbits, their migration reverses direction and they undergo orbital decay 
(Lin Sz Papaloizou 1985). The migration of some planets are stalled at intermediate radii when 



-31 - 



the disk is depleted. Note that the number of planets which stalled with intermediate semi major 
axes or cquivalently periods is a function of T(jcp/Tmig,i- The initial migration time scale at lAU, 
Tmig^i, provides an arbitrary fiducial comparison between the disk depletion and planets' migration 
time scales. In models with rdep/Tmig,i > 1, delicate timing is needed for a few planets to stall 
their orbital evolution at a fraction of an AU. This result is similar to those found by Trilling et 
al. (2002) and Armitage et al. (2002). But, in models with relatively small Tdep/rmig,!, a significant 
fraction of the gas giant planets formed near 1 AU may attain asymptotic semi major axes between 
0.04 to 1 AU. 

In Figure 11, we highlight the critical migration paths of (1) those planets which marginally 
migrate from an initial semi major axis ai to 0.04AU and (2) those planets which attain an asymp- 
totic semi major axis which are greater than 0.9 times that of their initial semi major axis 02. 
Since planets form with semi major axis between ai and 02 attain asymptotic a's between 0.04 AU 
and ~ 02, relatively large values of Tdep/Tnug,i would lead to a significant modification in the semi 
major axis distribution of the planets. 

5.3. The modification of the semi major axis distribution by migration 

We apply the above prescription for migration to the planetary growth model. We added the 
effects of migration to the Monte Carlo calculations in §4.2 (Figures 9 a, b, and c). We assume that 
type II migration starts when a planet's mass reaches Afg ^h and a gap is opened up in the nascent 
disk near its orbit (models (a) and (b)). In the case that the gas truncation is induced by isolation, 
it is less clear when type II migration would be initiated. For simplicity, we assume that migration 
starts when planetary mass reaches Mg^iso which is larger than Mg^th- For model (c), migration is 
initiated when tr becomes larger than h/1.5 such that there is a range for which gas accretion and 
migration occur concurrently. This migration condition is similar to Mp > Mg.vis- (As mentioned 
in the last section, the gas truncation by Mgas,vis seems to be inconsistent with the observational 
data, but the migration condition by Mgas,vis may be reasonable.) 

In these calculations, Tdcp = 10^ — 10^ years. The time dependent calculation of disk evolution 
(Lynden-Bell & Pringle 1974) indicates that the disk mass declines on the viscous diffusion time 
scale near R^a- If gas depiction in disks is due to their viscous evolution, we would expect Tdep to 
be comparable to rdisk,acc (eq. [70]) near ~ lOAU. In order to match the observed properties 
of protostellar disks around classical T Tauri stars, we adopt a = 10~^ which corresponds to 

Tdep/7disk,acc ~ 1 at lOAU. 

The results of our simulations are shown in Figs. 12 for three series of models. In each case, 
the gas and core accretion are truncated by the conditions which correspond to those in Fig. 9, 
respectively. The results show that the spatial distribution of the gas-poor cores is not affected 
by the migration because it only affects those planets which are able to accrete gas and to open 
up gaps. But for gas giant planets, equation (65) indicates that the migration time scale increases 
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with their masses and semi major axes. The less massive gas giants are formed preferentially with 

relatively small semi major axes and they migrate to ~ O.OAAU in all the cases. This result is 
consistent with the observed mass distribution of the short-period planets which appears to be 
smaller than that of planets with periods longer than a few months (Udry et al. 2003). 

Gas giant planets with Tmig ^ Tdisk migrate over extended radial distance provided the disk 
gas is preserved for a sufficiently long time for them to form. For example, the critical value of 
/disk for the formation of gas giants is ~ 3 — 8 at a ~ lAU where rjice = 1 (see §4.1). Prom 
eq. (18), we find that in disks with /disk larger than the critical value, the time scale for a core to 
acquire a mass Mp > 4M0 is 3rc,acc i$ lO^yrs which is much smaller than the mean disk depletion 
time scale Tdisk ~ 10^ ~ 10^ years. If the subsequent gas accretion is limited by either Mg or 
-^^g.iso) the asymptotic mass of the gas giant would be < lOMj (Figs. 9). For this range of Mp 
and /disk, we find, from eqs. (65) and (70), that the migration time scale Tmig is less than a few 
Myr provided a ~ 10~^. In most disks, Rm 3> lAU so that the planets formed at o ~ lAU have 
TcjOCC ^ ^mig ^ ''"disk- Thus, within our model here, most gas giant planets formed within ~ l^C/ 
are removed by the orbital migration independent of the truncation condition for gas accretion. 

Nevertheless, there is a population of gas giants which attain asymptotic semi major axis in 
the range of a fraction of an AU. In fact, the distribution of semi major axis is nearly continuous 
from ~ 0.04 — 4AU. As we have indicated in the discussion of Figure 11, planets formed between 
ai and 02 slightly outside the initial radius of maximum viscous stress (Rm) migrate outward first. 
They reverse their direction and migrate inward after iim has expanded beyond their orbit. (In all 
our calculations, we adopt R^ = lOAU initially so that ai ~ 3— 5AU and 02 ~ lOAU.) Their period 
distribution in the range of a few weeks to years is determined by the ratio of the disk depletion 
and viscous time scales. With our adopted value T(jep/Tdisk,acc(— 1); the simulated a distribution is 
similar to that observed. 

Numerical simulations of disk-planet interaction suggest that there may be lingering gas ac- 
cretion after the gap formation which may lead to Mp to be somewhat but not substantially larger 
than Mg^th- The post formation accretion may be particularly effective for the eccentric planets 
because their large epicyclic excursion can be extended beyond the boundaries of the gap. From 
eq. (65) we find that Tmig — ^disk when Mp reaches 

If such an effect leads to the formation of any massive planets with Mp > Mmig, they would not be 
able to migrate significantly from their birth place. 

5.4. Preservation of the planet-deficit domain 

The results in Figures 12 show that the domain of the planet desert becomes more clearly 
defined with all three prescriptions for the truncation of gas accretion. In comparison with the 
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planets' initial properties (see Figs. 9), the effect of the orbital migration is to reduce and erase 
the differences introduced to the Mp — a distribution by the different truncation conditions for gas 
accretion. 

In this series of simulations, we only consider type II migration, which occurs only if a planet 
can open up a gap with a mass in excess of (0.3-l)Mg^th [~ (30 — 10^) x (a/lAU)'^/^M®]. The 
lower- Mp boundary of the planet desert is not modified by the type II migration process because 
it is populated by cores which have not yet initiated the process of rapid gas accretion and may 
not migrate. As we have already shown that this boundary is primarily determined by the cores' 
mass required for transition from the predominantly planetesimal to the mostly gas accretion, it is 
insensitive to the truncation condition for the cores and their asymptotic masses. 

A direct comparison between Figures 9 and 12 clearly indicates that the planet migration 

process significantly alter the upper- Mp boundary of the planet-deficit domain. Gas giants formed 
with a < ai 3 — 5AU (sec Figs. 12) migrate to 0.04 AU. But some planets formed between ai and 
a2 ~ lOAU, attain asymptotic a between 0.04 — 4AU. As the semi major axis of these marginally 
survivable gas giants spread out, they attain a nearly logarithmic period distribution similar to 
the observed period distribution. Outside a2, the gas giants form in the advanced stages of disk 
evolution when there is insufficient gas to induce any significant migration. Based on the results in 
Figures 12, we infer a substantial upturn in the semi major axis distribution of Jupiter-mass gas 
giants at around lOAU. The boundary of the upturn in the period distribution is determined by 
the ratio of gas depletion and the viscous evolution time scales of the disk. 

During the migration, the gas giants are not expected to lose mass. Thus, the upper-Mp 
boundary of the planet-deficit region is set by the lowest gas-accretion-truncation mass between ai 
and 02- Above this upper boundary, the spread in Mp among the gas giants with any given a is 
due to the dispersion in the gas-accretion-truncation condition induced by the mass range of the 
host stars and nascent disks. 

In this series of simulations, the outer-a boundary of the planet desert in the Mp — a diagram 
occurs at ~ 3 AU for all three classes of models. The magnitude of the outer-a boundary is 

determined by the location where cores can attain a sufficient mass to undergo efficient gas accretion 
prior to the depletion of the gas in the disk. Since this condition is associated with the evolution 
of the cores, the location of the outer-a boundary does not depend on either the gas-accretion- 
truncation conditions or the migration process. 

Beyond the outer-a boundary, ice-giant planets, similar to Uranus and Neptune, can form with 
10 — IOOMq. In these region, the cores grow relatively slowly. Even though the asymptotic mass of 
the cores is relatively large, they are unable to acquire Mc^kh ~ 4A/0 to initiate a phase of rapid 
gas accretion. The large asymptotic mass of these planets is acquired during the final phase of 
their assemblage after gas in the disk has been depleted. In the limit that some residual gas in the 
disk is preserved for an extended length of time, type II migration cannot relocate these planets 
because their mass is too small to open up a gap. Nevertheless, the mass spectrum at any a beyond 
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the outer boundary of the planet desert can evolve slightly with gas accretion and significantly 

via residual core accretion and giant impacts. The time scale of growth through giant impacts 
is similar to that for the oligarchic cores and is an increasing function of the planets' mass. This 
further coagulation process leads to a continuous mass spectrum. The range of ice giants' masses in 
the region just beyond the outer-a boundary is determined by the dispersion in the heavy element 
content, i.e. /disk, out there. 

The center of the planet-deficit domain is not entirely void of planets with intermediate Mp 
and a's. Within the present theoretical framework of our models, only type II migration can modify 
the planets' semi major axis. But the necessary condition for type II migration is for gas giant to 
reach their gas-accretion-truncation mass -Mg th which is an increasing function of a. Gas giants 
with Mp less than the local Mg^th can be formed interior to any disk location. However, planets 
formed interior to the radius of maximum viscous stress in the disk must always migrate inward. 
Since i?m is generally comparable to or larger than 02 , planets with intermediate Mp and a cannot 
migrate into the planet-desert domains from smaller or larger a's. In our simulation, a few planets 
located in the planet-deficit, intermediate Mp and a region can be formed in marginally massive 
disks where the /disk > /disk.min a-nd /g ~ (Mp/Mj)(a/1AU)~^/^ (see eqs. [55] and [47]). In this 
estimate, we assume that since these planets do not have adequate mass to open a gap, they 
eventually acquire all the nearby residual disk gas. 



5.5. The spatial extent of planetary systems 

In the scenario we presented here, we artificially halt the migration at 0.04 AU to represent 
the effect of tidal and magnetospheric interaction between short-period planets and their host 
stars. Besides these processes, MHD turbulence in the disk may also play a role in halting their 
migration (Terquem 2003). Since most of the short-period planets were formed with a < ai, they 
generally form with less mass than the planets which accreted gas later, at larger disk radii, and 
retained a greater fraction of their original semi major axis. Most gas giant planets formed interior 
to ai migrate to 0.04 AU. Since gas giants can readily form well within ai (see Figs. 9), a large 
accumulation of short-period planets is expected. Although there is a spike in the observed period 
distribution just outside 0.03AU, the amplitude of this enhancement is modest. Near the surface 
of the host stars, the combined effects of tidal, magnetospheric, and radiative effects of the host 
stars as well as the secular perturbation of more distant planets, disks, and stars may lead to the 
destruction of the short-period planets and the disruptions of their orbits. 

For the outer boundary of planetary systems, we find that there are several very massive 
planets which managed to migrate to a ~ 50 — lOOAU in panels a), b), and c) of Figure 12. 
These planets were formed at 10-30 AU, just outside 02 in massive disks we have simulated. With 
their large values of a ~ 10-30AU and /disk ~ 20-30, the gravitational stability parameter Q 
(eq. [56]) at formation locations of these planets are found to be 0.7-1.2, which may be (marginally) 
gravitationally unstable. Hence, the outward migrating massive gas giants may actually be more 
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rare. 

For disks with modest masses, ice giants can accrete gas to form sizable gaseous envelopes 
(represented by open blue circles) well beyond the ice boundary. Although we assume that the ice 
giants do not migrate because they do not have a sufficient mass to open a clean gap. Nevertheless, 
the weak tidal interaction may be adequate to induce the ice giants to migrate primarily outward 
over a limited extent. In the Solar System, the migration of Neptune over several AU's on a time 
scale of ~ 10 Myr may have led to the capture of Kuiper Belt objects on its 3:2 (but not on its 2:1) 
mean motion resonances (Ida et al. 2000). 

In the current analysis, we do not consider the interaction between multiple planets which can 
lead to dynamical instabilities (Chambers et al. 1996; Lin k, Ida 1997; Marzari & Weidenschilling 
2002). Although sizeable planets may be ejected to large distances from their host stars, such 
an instability is likely to induce the orbits of the residual planets to become highly eccentric and 
dynamically segregated with large (> 10) period ratios, e.g., (Weidenschilling & Marzari 1996; 
Lin &L Ida 1997). The search for the orderly-packed ice giants with nearly circular orbits may 
provide useful constraints to differentiate the core accretion and gravitational instability scenarios 
for planet formation. 

6. Summary and Discussions 

We have investigated the mass and semi major distributions of extra solar planets through 
theoretical models based on the core accretion model for formation of gas giants. Core accretion 

from planetesimals is modeled by the results of N-body simulations. Gas accretion onto the cores is 
modeled with Kelvin-Helmholtz contraction of gas envelope. Since the termination of gas accretion 
onto cores have some uncertainty, we considered various models for the asymptotic mass of the gas 
giants. In some calculations, we also include the effects of type II migration. 

In general, the formation of gas giants is favored just outside the ice boundary at ~ 3AU. 
Inside the ice boundary, because the volatile gases cannot condense into grains, the mass of solid 
cores tend to be smaller than AM^ which is necessary for the onset of rapid gas accretion onto 
them. On the other hand, in outer regions (> lOAU), core accretion is so slow that cores cannot 
attain 4Mg before the disk gas is depleted on the time scale of ~ 10^ — 10^ years. In very massive 
disks, however, cores become massive enough for the onset of rapid gas accretion even inside the ice 
boundary within the disk life time to attain gas mass > lOOM®. Since planets grow very rapidly 
from IOM0 to IOOM0, planets of the intermediate masses are rare at a < 3AU ("planet desert"). 
If the effect of type II migration is included, the deficit becomes more clear and tends to become 
independent of the model of termination of gas accretion. It also results in accumulation of planets 
at a ~ 0.03 - O.IAU. As a result, the deficit becomes ~ 10 - lOOM® and 0.2 - 3AU. 

The current radial velocity survey of nearby stars is limited by the availability of accurate 
radial velocity measurements (down to the level of Vr > 10 ms"-*^) over the past few years. Planets 
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with masses comparable to that of Jupiter and a < 4-5AU are inferred from these observations. 

The processed data show a hint of the planet desert (Figure 1). Further observations in the regime 
of longer periods (larger a) and smaller Vr will clarify the existence of a planet desert (a deficit of 
planets with intermediate mass and semi major axis planets). 

Gas giants could also be formed by self-gravitational collapse of a disk, e.g., (Boss 2001). If 
extra solar planets are formed from the self-gravitational collapse, the mass and semi major axis 
distributions of extra solar planets should show very different features from those we show here based 
on the core accretion model for gas giants. As we showed in §4 and 5, with the core accretion model, 
gas giants may be limited to the regions a < 20AU, except for planets on eccentric orbits that are 
scattered from inner regions by dynamical instabilities in multiple gas giant systems (which we do 
not consider in this paper). On the other hand, in the self-gravitational collapse model, gas giants 
are preferentially formed far from their host stars (see eq. [56]), although the mass distribution of 
gas giants is not clear. The search for gas giants on nearly circular orbits beyond 20-30AU may 
provide useful constraints to differentiate the core accretion and gravitational instability models for 
formation of gas giants. 

If the the presence of a planet desert is confirmed by future observations, the boundaries 
of the planet-deficit domain in the mass and semi major axis plane should provide the following 
constraints on the planetary formation processes (see also Fig. 13): 

• The lower mass boundary (Mi) indicates the core mass that can initiate the onset of rapid 
gas accretion. Our theoretical model predicts that Mi ~ 4M0 and it does not depend on a. 
(If we include the coalescence between cores after gas depletion, Mi slightly increases with a.) 

• The value of upper mass boundary (M^) and its dependence on a provide clues on the 
truncation mechanism of gas accretion. For various models we find (1) Mu = 3 X 10^ - lO^Me 
(and independent of a) if the gas giants' growth is limited by an isolation process, (2) = 
2O(Q;/lO~^)(a/lAU)^/^M0 if a complete gap is formed through the viscous condition, or (3) 
Mu = 8O(a/lAU)^/^M0 if it is due to the thermal condition. However, if the effect of type II 
migration is included, the results in models (1) and (3) lead to similar Mp — a distribution. 
Then the magnitude of Mu constrains the migration timescales and disk depletion timescales. 

• The smaller a boundary (ai) indicates the maximum radius at which type II migration can 
be halted. 

• The larger a boundary (ou) traces out the radius where a dominant condition for formation 

of gas giants changes between rc,acc condition and Mc^iso condition. At a> a^, the former is 
more important than the latter, and vice versa. Our theoretical model predicts that Ou ~ 3 
AU. 

The results presented here show that the sequences of rocky (terrestrial), gas giant and ice 
giant planets are clearly segregated in the Mp — a diagram. The importance of this diagram to 
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planet formation is similar to that of the color-magnitude diagram for stellar evolution. Comparison 
of future observations and more refined theoretical models will reveal the formation processes of 
extra solar planets. In the near future our inability to directly detect the presence of terrestrial 
planets limits the comparison to only the giant planet population. In follow-up papers, we will 
consider a wider variety of models and identify some robust constraints. Our eventual objective is 
extrapolate the ubiquity of rocky planets and formulate a predictable theory of planet formation 
from the statical properties of the giant planets. 

We thank E. Kokubo, M. Ikoma, and G. Laughlin for useful discussions. This work is sup- 
ported by the NASA through NAGS5-11779 under its Origins program, JPL 1228184 under its SIM 
program, and NSF through AST-9987417. 
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Fig. 1. — Distributions of semimajor axes (a) and masses (Mp) of discovered extrasolar planets 
are shown by open circles. Unit of mass is Earth mass (Jupiter mass is Mj ~ 32OM0). % is 
the angle between an orbital plane and the line of sight. The dashed ascending lines correspond to 
radial velocity amplitude of 100 ms~^, 10 ms~^, and 1 ms~^ (assuming the host star mass is IM©). 



-43- 




Fig. 2. — as a function of a that we used, /d is an enhancement factor from the minimum-mass 
disk model for our Solar system. The jumps at 2.7AU are caused by ice grain condensation at 
a > 2.7AU. 
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Fig. 3. — (a) Analytical estimate of the core mass after 10^ years (upper panels) and 10^ years 
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Fig. 4. — The time evolution of core mass (10^, 10^, and 10^ years) is shown by sold lines, in the 
case of (a) /d = 1 and (b) /d = 10. Core mass is truncated by -/Vfc,iso (upper panels) or by -/Wc,noiso 
(lower panels), which are shown by the broken lines. The jumps at 2.7AU are caused by increase 
in solid materials due to ice condensation. 



-46- 





Fig. 5. — Core masses after 10^ years (upper panels) and 10'' years (lower panels) as a function of 
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Fig. 6. — Planet masses including gas envelope at 10^ years as a function of friisk (= fd = fe n) and 
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Fig. 7.— 



The the same plots as Figures 6 except for core truncation by Mr 
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Fig. 8. — The /disk distribution we use, which is a gaussian distribution in terms of log^o /disk with 
a center at log^o /disk = 0.25 and dispersion of 1. We omit the high /disk tail at /disk > 30, since 
such heavy disks are self gravitationally unstable. 



- 50 - 



6 



10* yiiiij — I 1 1 1 11 




rrjFFg 10* y i — i i iiiiii| — i 1 1 iiiii| — i mhb 



1 




ft [AU] G [AtJ] 

ICf* Niiii| — I 1 1 1 iiii| — I 1 1 iiiii| — I iji-ng 10^ N . — I I iiiiii| — I 1 1 iiiii| — I [jmH 





1 10^ Iff 
a [AU] 



a [AU] 



Fig. 9. — Theoretically predicted distribution based on core accretion model for gas giant planets 
(for the range of parameters we used, see text). Cores are truncated by Mc,iso- Gas accretion is 
truncated by (a) Mg^iso, (b) and (c) Mg th, and (d) Mg^vis- We adopt Aog = 2rH in (a), the critical 
Hill's radius rH,c being h and 1.5h in (b) and (c), respectively, and a = 10~^ in (d). The green 
filled circles and the blue crosses represent rocky and icy planets with gaseous envelopes less massive 
than their cores. The green and blue open circles represent gas-rich rocky and icy planets with 
gaseous envelopes which are one to ten times more massive than their cores. The red filled circles 
represent gas giants with envelopes more massive than ten times of their cores. For comparison, we 
also plot observational data of extrasolar planets in (d). The dashed ascending lines correspond to 
radial velocity amplitude of 100 ms^^ (upper), 10 ms~^ (middle), and 1 ms~^ (lower), respectively, 
assuming the host star mass is IMq. 
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Fig. 10. — The limits of masses of rocky and icy planets. The masses arc limited by the minimum of 
-^c.noiso (eq. [23]), Mc^sca (cQ- [27]), and Mc,iifc with t = IG years (eq. [60]). The solid, long-dashed, 
short-dashed, and dotted lines express the cases of /disk = 30, 10,3, and 1, respectively. 
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Fig. 11. — Planetary migration obtained by integrating dp with eqs. (65) and (68) is shown in sohd 
lines. Inside the thick dashed lines, final semi major axes reach 0.04AU. They are > 0.9 times their 
original locations, outside the thick solid lines. Dashed lines express Rm. Time is scaled by Tmig at 

lAU (rinig,l). (a) Tdep/Tmig = 10, (b) Tdep/Tmig = 1, (c) Tdep/Tmig = 0.1, and (d) Tdep/Tmig = 0.01. 
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Fig. 12. — Similar plots as Figures 9, but the effect of type II migration is included. The value of 
a-viscosity is taken as a = 10~^ to be consistent with disk depletion times ~ 10^ — 10'' years, (a) 
Gas accretion is truncated by -Mg^iso and core accretion by Mc^iso, (b) Mg^iso and Mc^noiso and (c) 
Mg^th and Mc^so- We adopt Aag = 2rH in (a), M^, = IMq in (c). 
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Fig. 13. — An illustrative figure for the planet desert. The lower mass boundary (Mi) indicates the 
core mass that can initiate the onset of rapid gas accretion. The upper mass boundary (Mu) and 
its dependence on a provide clues on the truncation mechanism of gas accretion. The smaller a 
boundary (ai) indicates the maximum radius at which type II migration can be halted. The larger a 
boundary (cu) traces out the radius where a dominant condition for formation of gas giants changes 
between Tc^acc condition and Mc^iso condition. For details, see text. 
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Mc 


a core mass 






Me 


an embryo mass 


see 


§2.4 


Mc,iso 


a core isolation mass 


eq. 


(20) 


-Mc^noiso 


a maximum mass of a core 
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(22) 


Me,iso 


an embryo isolation mass 


eq. 


(25) 


Me^sca 


an embryo mass limited by scattering 


eq. 


(27) 


Me,iife 


an embryo mass at t = r* (e,g., IG years) 


eq. 


(60) 


Mc^crit 


a critical core mass for initiating gas accretion onto the core 


eq. 


(33) 


Mc,KH 


a core mass for K-H contraction time being = rjep 
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(31) 


Mp 


a planet mass including core and envelope 






Mg^iso 


a gas giant isolation mass 
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(35) 


Mg^noiso 


a maximum mass of a gas giant 


eq. 


(47) 


^g,th 


a gas giant truncation mass by thermal condition 


eq. 


(38) 


Mg,vis 


a gas giant truncation mass by viscous condition 


eq. 


(37) 


Mg,feed 


a supplied gas mass into a feeding zone by viscous diffusion 


eq. 


(46) 


''"c,acc 


a core accretion time scale before gas depletion 


eq. 


(6) 


'''e,acc 


an embryo accretion time scale after gas depletion 


eq. 


(9) 


'''dep 


a depletion time scale of disk gas 


eq. 


(48) 


TKH 


a Kelvin-Helmholtz contraction time scale of gas envelope 


eq. 


(34) 


T'mig 


a migration time scale of a giant planet 


eq. 


(65) 


Sd 


dust surface density of a disk 








gas surface density of a disk 






^g,o 


initial gas surface density of a disk 


eq. 


(48) 




an enhancement factor of Sd due to ice condensation 


eq. 


(11) 


/d 


an enhancement factor of Ed from a nominal disk 


eq. 


(11) 


/g 


an enhancement factor of Eg from a nominal disk 


eq. 


(12) 


/g,o 


an initial enhancement factor of Eg from a nominal disk 


eq. 


(49) 


/disk 


= /g,o = fd in this paper 






/tg 


a value of /disk to segregate terrestrial/gas-giant planets regions 


eq. 


(50) or (51) 


/gi 


that to segregate gas/ice giants regions 


eq. 


(52) 




ice boundary (a semi major axis of ice condensation) 


eq. 


(15) 


«tg 


a semi major axis to segregate terrestrial/gas-giant planets regions 


eq. 


(53) 


«gi 


that to segregate gas/ice giants regions 


eq. 


(54) 



Table 1: List of notations 



